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Abstract 

Electromagnetic form factors of the nucleon from relativistic quark models are 
analyzed: results from null-plane projection of the Feynman triangle diagram 
are compared with a Bakamjian-Thomas model. The magnetic form factors 
of the models differ by about 15% at spacelike momentum transfer 0.5 GeV 2 , 
while the charge form factors are much closer. Spurious contributions to 
electromagnetic form factors due to violations of rotational symmetry are 
eliminated for both models. One method changes magnetic form factors by 
about 10%, whereas the charge form factors stay nearly the same. Another 
one changes the charge form factor of the Bakamjian-Thomas model by more 
than 50%. 
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I. INTRODUCTION 



In light-cone quark models wave functions may be boosted kinematically, that is, inde- 
pendent of interactions. Because form factors involve boosted wave functions, this attractive 
feature has motivated many electroweak form factor calculations |],|2| in Dirac's light-front 
form of relativistic few-body physics 0]. When the quark model is Lorentz invariant and 
the exact (interaction dependent) current is evaluated, the form factors will be Lorentz 
invariant. Because most form factor models are based on free quark currents (in impulse ap- 
proximation, IA, and/or other approximations), the resulting form factors are in fact frame 
dependent. 

The light-front form is obtained from the instant form in the infinite momentum limit, 
and this amounts to the well-known change of momentum variables (p + = p° +p z ,p~ = 
p° — p z ) on the light cone [[|. Here the conventional choice of the light cone direction is 
the z-axis and the light cone time is x + = t + z. Because in front form both transverse 
rotation generators become interaction dependent, rotational symmetry is more difficult to 
implement. Thus light front models include kinematic boosts at the expense of interaction 
dependent angular momentum operators. 

Light front models obtained from quantum field theories usually involve interacting spin 
operators. There is a unitary (but interaction dependent) transformation to a representation 
where spins are free and state vectors are specified on a null plane u • x — with a lightlike 
(i.e. uj 2 = 0) four- vector uj^ 1 ||. If a theory is Lorentz invariant and the exact current 
operator is known, e.g. in quantum field theory given by the Mandelstam prescription, its 
form factors are expected to be independent of u; M . In this case the electromagnetic current 
operator in general has many-body contributions |J. When it is approximated by a sum of 
free one-body currents, as in the impulse approximation for example, Lorentz invariance of 
the form factors is violated. 

When Lorentz invariance of the form factors is violated, additional cu-dependent current 
components arise along with their form factors which are needed to fully parameterize the 
current matrix elements that now obey fewer constraints. When Lorentz invariance of the 
quark model form factors is fully restored by including two-body currents, these spurious 
currents are canceled and the remaining calculated form factors are modified as well. In 
other words, achieving cu-independence of form factors by eliminating spurious from factors 
is just one step towards Lorentz invariance but, by itself, it does not guarantee correct form 
factors yet. In the case of a»-dependence the unphysical currents can be removed. As a bonus, 
the inconsistencies between current matrix elements of different helicities are reduced. This 
is done here for electromagnetic nucleon form factors of constituent quark models in light 
front dynamics. These models are kept as simple as possible (e.g. no anomalous magnetic 
moments of constituent quarks) because we are not attempting to fit the data. Rather we 
wish to compare their electromagnetic form factors and the size of errors due to violations 
of the rotational symmetry. 

Two light-cone versions of the constituent quark model (NQM), one based on the 
Bakamjian-Thomas prescription pj, called BT in the following, and another on the null- 
plane projection of a triangle Feynman diagram (see Araujo et al. in f2|), called KW (see 
Konen et al. in || for a = 1/2), are briefly reviewed in Sect. II. The electromagnetic cur- 
rent matrix elements are given in Sect. Ill, and the results in Sect. IV and conclusions are 
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presented in Sect.V. 



II. LIGHT-CONE QUARK MODELS 

We start with a brief description of the light front notation. The longitudinal quark 
momentum fractions are defined / P + with the total nucleon momentum P + = 

J2iPt so that Yli x i — 1 an d < Xi < 1, where Pi are the quark momentum variables. The 
ith quark momentum in the nucleon rest frame is given by 

h± = Pi-L - XiP±, (1) 

— * — * . 

with k± = (k x , ky) etc. so that these internal variables satisfy J2i ki± = and kj~ = 0. For 
a three-quark bound state the relative four-momentum variables are the space-like Jacobi 
momentum variables in which the kinematic invariants x« play the role of masses. For 
example, is the relative quark momentum between the up quarks of the proton in the 
uds-basis and Q3 between the down quark and the up quark pair, so that for the + and 
_L= (x, y) components 

<?3 = ; , Qi = {xi + x 2 )P3 - :c 3 (pi + p 2 ), (2) 

x x + x 2 

etc. Because q% = = Q% both relative momentum variables are space-like. In the uds- 
basis II quarks are treated as distinguishable, and up and down quarks are symmetrized 
explicitly. The antisymmetric color wave function is suppressed. For the proton (neutron) 
the down (up) quark carries the label 3. 

In light front dynamics the total momentum motion rigorously separates from the internal 
motion. Therefore, the internal nucleon wave function ip(x{, q%, Q3, Xi) does not change under 
kinematic Lorentz transformations or translations. Here the Xi are the quark helicities. Thus, 
if the wave function is known in the nucleon rest frame it is known everywhere because the 
seven kinematic generators are transitive in the null plane. 

Three-quark wave functions for baryons have been constructed in light cone versions of 
the constituent quark model (LCQM) to study the static properties of the nucleon |[],[| and 
electromagnetic transition form factors for iV — * iV* and iV —>■ A processes P,[Hl|j. Such 
relativistic quark models significantly improve many predictions of the NQM of which the 
nucleon weak axial charge gA is the best known example. For a recent analysis of the static 
observables and nucleon electroweak form factors see ref . |L2| . 

The three-quark wave function for a nucleon is the product of a totally symmetric mo- 
mentum wave function and a nonstatic spin wave function xx which is an eigenfunction of 
the total angular momentum (squared) and its projection on the light cone axis. The spin 
wave function can be represented as a linear combination of products of matrix elements 
between valence quark light-cone spinors coupled by appropriate 7 (and isospin)-matrices 
to the spin and isospin of the nucleon. Free Melosh rotations JTTJ] are of central impor- 
tance in LCQMs for the construction of relativistic many-body spin-flavor wave functions 
for hadrons from those of the nonrelativistic quark model (NQM). Instant Dirac-spinors are 
transformed into light-front helicity eigenstates (light-cone Dirac-spinors) by Melosh trans- 
formations which include the kinematic quark boosts. Direct coupling of Pauli spinors by 
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SU(2) Clebsch-Gordan coefficients leads to a set of spin-flavor wave functions that are in 
one-to-one correspondence with the NQM states. (For details, see ref. [0.) 
The wave function of the nucleon with nucleon helicity A is here taken to be 

ip N (\) = N<j)(x 3 , q 3 , Q 3 )Xa, (3) 

where N is a normalization constant fixed by the proton charge (electric form factor at q 2 = 
0). The Gaussian momentum wave function is written in terms of the relative momentum 
variables according to the Brodsky-Huang-Lepage prescription WA 



x 



q3,Qs)=e~ M ^ 2 (4) 



with a size parameter f3 and the free mass operator M 3 , i.e. the sum of the quark light-cone 
energies in the nucleon rest frame, 



mi = y *£- + m ? = - q \ 1 —^ - Q " - +y^. (5) 

" , X Xi x x x 2 x 3 (l - X 3 ) ^ Xi 



In the projection of the triangle Feynman diagram to the null plane the totally symmetric 
nonstatic proton spin wave function \\ ( see [0 f° r more details) originates from the quark- 



nucleon Lagrangian |]12 | 



£ q *~N = (qi[id N } l5 iT 2 Cql) ■ q 3 V N + (23)1 + (31)2, (6) 
where and ^jy are the quark and nucleon fields, C = ry 2 7° is the charge-conjugation 



operator, and 75 is characteristic of the spin of the quark pair ]T6| . The term [P] = 1 '^^ N 
with P the total momentum and the nucleon mass for the null-plane projection model 
(KW), while for the Bakamjian-Thomas model (BT) the nucleon mass in [P] is replaced 
by the free three-quark mass M 3 of Eq. |^. In the BT-model [P] derives from the Melosh 
rotation of the quark Pauli spinors. The nucleon positive-energy projection operator [P] 
corresponds to the mixing parameter a = 1/2 in Araujo et al. of [0 of two of the three 
linearly independent three-quark-nucleon couplings that reduce to the NQM S-wave ground 
state in the static limit [|13|,[[^]. We restrict ourselves to the case a = 1/2 because it can be 
directly compared with the BT-model. 

In the triangle diagram for the form factors the quark fields are replaced by light-cone 
spinors Ui, quark 3 being the down quark of the proton in the uds-basis, so that the con- 
ventional nonstatic spin wave function 

Xn = u 2 [P]^Cul ■ uiU N - u 3 [P]^Cul ■ u 2 u N (7) 

results from Eq. []. The momentum wave function (M| — m^)0 in Eqs. |3],|] is proportional 
to the vertex function of C q z_ N of Eq. || 

These simple relativistic quark models depend on two parameters, the common u, d 
constituent quark mass m q and the size constant f3 of the confinement potential. The 
proton mean square radius determines f3 so that 1//3 ~ (r 2 )p 2 up to relativistic corrections. 
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III. ELECTROMAGNETIC FORM FACTORS OF THE NUCLEON 



If J M is the full electromagnetic current, the current matrix element of a nucleon state 
consists of the Dirac and Pauli conserved currents, 



(N(P>, X')\J»\N(P, A)) = u N (P', A')[ 7 ^i(<7 2 



2m 



N 



F 2 (q 2 )}u N (P,X), 



(8) 



where <r M1/ = z[7 M , 7"]/2. We will also use the standard electric and magnetic Sachs form 
factors 



G E {q 2 ) = F^q 2 ) - V F 2 (q 2 ), G M = F 1 + F, 



2- 



(9) 



where r] = —q 2 /Am 2 N > 0. 

If free quark currents are used as in the impulse approximation, the possible in- 
dependence of the electromagnetic current matrix elements between baryon states leads 
to three additional conserved components [17] with form factors so that Eq. ||] becomes 

F-2 



(N(P',X')\J2e j q J l%\N(P,X))=u N (P',X') (F l7 " • i-^v' "v(P.A) 



+u N (P', X')B X I ' : 



1 



(O • P (1 + i])m N , 
+u n (P',X')(b 2 ^-^ + B 3: ^ 



(10) 

(P^ + P^)u N (P,X) 



■uo-P 



;U • ^ U N (P,X), 



where u ■ q = q + = for the standard choice of light cone axis. 

In front form all form factors are usually extracted from matrix elements of J + in the 
light-cone analog of the Breit frame (where q + = 0), the so-called 'good' current compo- 



nent 18 



eq 



2m 



eF^q 2 
-F 2 (q 2 



m N 
P+ 
m N 



A' 



P 



(NiP'W+lNiP),}, 



(11) 



with q L = q x — iq y . This is no longer possible when cu-dependent currents are present |L9 
We therefore define the more general current matrix elements 

J^ X = (N',X'\J^N,X) (12) 

in the rest frame of the nucleon. Just as Eq. [11] follows from Eq. |8| we obtain 

1 

a 





= P" 




+ 2B 1 




\m N 




]+ 


= P" 


•'(- 


F 2 
2m 2 N 


T L 


= q L 


c 


f F2 


\mN 


2mN 



m N 
aBi 



+ B l 



— a 



m N 



jL 

J n 


= {q > 


J i1 


■A 







Bid \ 



Am 2 N 2rriN J 



F 2 B x a 



4m 2 N ' 2m N 



Bi 



2m 



N 



m N 



(13) 



from Eq. |10| and [13|, where a = 1/(1 + r])m^. It is straightforward to check the symmetries 

(14) 



J ii - TT' 



L ' 



etc. From Eq. [T^ in conjunction with Eq. [13| we find two solutions for the Dirac and Pauli 
form factors 



j+ 

Fi = m N {(2am N - 1)^1 + 2m 2 N 



— a 



-m%{ 



2a J 



TT 



+ 2 



m N P + [m N 
which we label set 1, while 



— a 



m N 
X 



P+q 

2d J || 
p+qL rriN q L 



J+ jL 

+ 2(1 -am N )^}, 



(f 



Tl 



}■ 



(15) 



— a 



-J 



j+ 

Fx = m N {^ + 4m 
l 2 = -i m %{[l-am N }^ + ^}. 



JL 

-l-l j Ti- 



m N 
J_ L 

iJ L ) 2 ' 2 q L 



\2m N q L (q L ) 2 



(16) 



is the set 2. 

Subtracting the spurious form factors according to the prescriptions of Eqs. |15|,|16] will 
be called corrected form factors in the following. 



IV. NUMERICAL RESULTS 

The form factor calculations involve matrix elements of products of up to ten Dirac 
spinors which are evaluated and summed using symbolic codes generated by Mathemat- 
ica. For the KW-model the transverse momentum integrations are carried out analytically, 
so that only the integrations over the longitudinal momentum fractions x\,X2 are evalu- 
ated numerically. The BT-model involves six-dimensional numerical integrations which are 
performed by Monte Carlo routines |[20|| . 

Although relativistic models yield predictions for any value of the momentum transfer, 
the constituent quark concept is expected to lose validity (i.e. the momentum dependence of 
the quark selfenergy becomes less important) at the spontaneous chiral symmetry breaking 
scale (x = 47r/ w ~ lGeV), where current quarks are expected to become the relevant degrees 
of freedom. This limits the validity of the models to fairly small values of Q 2 = —q 2 . Both 
of our light-cone quark models are based on pointlike constituent quarks. (Better fits can 
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be obtained using small anomalous magnetic moments and a finite size form factor, but this 
is not our objective here. See Cardarelli et al. in [|[].) 

In figures 1 to 4 we compare the electromagnetic form factors of the KW and BT models 
(uncorrected and corrected) for the same parameters (3 = 0.32 GeV and m q = 0.32 GeV. 
Results are shown relative to the dipole form factor Fp = (1 — g 2 /0.71 GeV 2 )~ 2 with q 2 in 
GeV 2 except for Fig. 3 (neutron charge form factor). While the charge form factors of KW 
and BT (both uncorrected) for the proton in Fig. 1 are close to each other below 0.7 GeV 2 , 
both models' magnetic form factors in Fig. 2 differ by a roughly constant amount 0.2, i.e. 
20% at 1 GeV 2 . The corrected magnetic proton and neutron form factors BT1 and BT2 
in Figs. 2 and 4 are close to each other but differ by about 20% from the uncorrected BT- 
model. This we interpret to mean that the corrected magnetic form factor of the BT-model 
is nearly unambiguous and correct. The corrected KW1 and KW2 form factors differ by a 
few percent (Figs. 2 and 4), however the correction itself is smaller than in the BT case. 
The charge form factors KW1 and KW coincide, so that the dotted line is not visible in the 
figures and KW2 leads to a slight correction in the proton case and a larger correction for 
the neutron. For the proton charge form factor BT1 and BT2 differ much more, so that it 
is poorly predicted by the BT-model in contrast to the KW-model. 

Better fits can be obtained for both models, but for different parameters sets, and samples 
are shown in Figs. 5-9. In contrast to the NQM, the proton charge form factor G P E falls faster 
than G P M in both relativistic quark models. This is a general feature of LCQMs, which is in 
qualitative agreement with the recent data from Jefferson Laboratory ||22|| . Yet the proton 
charge form factors of both models fall off too fast compared to the JLab data in Fig. 5. 

The major disagreement involves the neutron charge form factor that is low by a factor 
of more than 2. However, this is not surprising because it is known to be sensitive to the 
spin force and pion cloud effects that are beyond the scope of the simple relativistic quark 
models considered here. 

The trends of uncorrected form factors persist when spurious form factors are subtracted 
according to Eqs. p~5| , p76] as shown in Figs. 1 to 9. From Fig. 1 we see that the proton charge 
form factors of the KW-model and the BT model for set 1 are hardly affected by the spurious 
contributions, while the corrections for the BT-model for set 2 are more than 50%. For the 
magnetic nucleon form factors in Figs. 2, 4 the corrections are similar for both models. 



V. SUMMARY AND CONCLUSION 

The impulse approximation leads to violations of gauge invariance and rotational sym- 
metry by LCQMs, i.e. spurious dependence on the choice of light-cone axis. Depending on 
which helicity matrix elements are used to extract electromagnetic form factors, there are 
different schemes of subtracting spurious components from these form factors. Subtracting 
spuriosities reduces inconsistencies between helicity matrix elements. When electromagnetic 
form factors roughly agree that are obtained from different correction schemes, then we con- 
sider them as essentially correct and reliably predicted by the LCQM in question. This 
is the case for the magnetic proton and neutron form factors of the KW- and BT-models. 
The KW-model also predicts G P E reliably in this sense, in contrast to the BT-model that 
predicts G P E poorly. The typical corrections for spuriosities are larger for the BT-model as a 
consequence of its larger high-momentum components that enter via the three-quark mass 
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M% in the quark-nucleon coupling in Eq. |6|. (Note that M3 of Eq. |5| depends on the internal 
quark-momentum variables.) For — q 2 < 1 GeV 2 , the corrected form factors differ from the 
uncorrected ones by about 20% for the KW-model and up to 50% for the BT-model. These 
are also the typical differences between reliably predicted form factors of these models. Thus, 
the theoretical errors of the form factors of the models are of roughly the same size as the 
differences between reliably predicted form factors of the models, so that it is difficult to 
distinguish them. These findings suggest that it is important for LCQMs to improve the 
impulse approximation and use conserved electromagnetic currents that are consistent with 
the interactions of the models. 
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FIGURES 
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FIG. 1. Proton charge form factor G P E /Fn relative to dipole shape Fr> for m q = (3 = 0.32 GeV: 
The solid line is KW uncorrected while dashed and dotted lines KW1, KW2 are the corrected set 
1,2 cases, respectively. The BT-model curves are similarly denoted by BT, BT1, BT2. 
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FIG. 2. Proton magnetic form factor G P M / HpFu relative to dipole shape fJLpFjj. The notation 
and parameters are as in Fig. 1. 
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FIG. 3. Neutron charge form factor with the notation and parameters of Fig. 1. 
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FIG. 4. Neutron magnetic form factor / fi n F£> relative to dipole shape h p Fd with the no- 
tation and parameters of Fig. 1. 



14 




FIG. 5. Fit of G P E /F D for the KW-model (KW) with m q = 0.34 GeV,/? = 0.32 GeV and 
BT-model (BT) with m q = 0.26 GeV,/3 = 0.36 GeV in the notation of Fig. 1. The results 
corrected for spurious form factors according to sets 1,2 are denoted by KW1, KW2 and BT1, 
BT2. The full circle data are from ref. [21], the full squares are from [22]. 
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